Collective states in nuclei and many-body random interactions 



R. Bijker^ and A. Frank^^^ 
^ Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 
^ Centro de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 139-B, Cuernavaca, Morelos, Mexico 
(April 3, 2000) 

Low-lying collective states in nuclei are investigated in the framework of the interacting boson 
model using an ensemble of random many-body interactions. It is shown that whenever the number 
of bosons is sufflciently large compared to the rank of the interactions, the spectral properties are 
characterized by a dominance of — 0^ ground states and the occurrence of both vibrational 
and rotational band structures. This indicates that these features represent a general and robust 
property of the collective model space. 
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I. INTRODUCTION 



Random matrix ensembles provide a powerful tool to study generic spectral properties of complex many-body 
systems Most applications in the literature have centered on global characteristics such as first neighbor energy 
distributions, which typically involve states with the same quantum numbers (angular momentum, parity, isospin, 
...). Recently, the relation between low- lying states in even-even nuclei with different quantum numbers was examined 
using Hamiltonians with random interactions in the nuclear shell model (SM) |^-|^ and the interacting boson model 
(IBM) 1^,0. These studies have given rise to several surprising results. In both cases it was found that for a large 

variety of conditions there is a dominance (~ 60%) of = 0+ ground states despite the random nature of the 
interactions. In addition, in the SM strong evidence was found for the occurrence of pairing properties Q, and in the 
IBM for both vibrational and rotational band structures These results are not only based on energies, but also 
involve the behavior of the wave functions via the pair transfer amplitudes in the case of pairing, and the quadrupole 
transitions for the collective bands. The use of random interactions (both in size and in sign) show that these regular 
features arise for a much wider class of Hamiltonians than are usually considered to be realisitic. These results are 
in qualitative agreement with the empirical observations of very robust features in the low-lying spectra of medium 
and heavy even-even nuclei and a tripartite classification in terms of a seniority, a vibrator and a rotor regime [^,^ . 
The conventional wisdom in nuclear structure physics is that the observed properties of nuclei can be explained by 
specific features of the SM (or IBM) Hamiltonian. The studies with random interactions, however, seem to imply 
that some of the generic characteristics of these systems may already be encoded in the corresponding shell model (or 
sd boson model) space. This is particularly striking in the case of the IBM, for which the model space corresponds 
to a drastic truncation of the original (shell model) Hilbert space to that composed of like-nucleon pairs with angular 
momentum L = and L = 2 |10 1^. The selection of such a restricted subspace seems to impose strong constraints 
on the possible spectral properties. 

These considerations lead naturally to the question of what arc the specific causes of this behavior, given that 
the ingredients of the calculations arc the preservation of fundamental symmetries of the Hamiltonian (hermiticity, 
rotational invariance, time-reversal invariance), the one- and two-body nature of the Hamiltonian, a given number of 
active particles, and the structure of the model space. In it was shown that the preponderance of = 0+ ground 
states in the nuclear shell model is not due to the time-reversal symmetry of the interactions. The purpose of this 
paper is to address explicitly the role of the particle number and the rank of the random many-body interactions on 
the systematics of collective states in nuclei. 



II. RANDOM INTERACTIONS IN THE IBM 



To study the global features of low-lying collective states in nuclei we carry out an analysis of the IBM with random 
interactions. In the IBM, collective nuclei are described as a system of N interacting monopole and quadrupole bosons 
. We consider all possible one-, two- and three-body interactions. The one-body Hamiltonian contains the boson 
energies 
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The two-body interactions can be expressed as 
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with Plm = (— 1)^ ^-Pl,-m- Here denotes the creation operator of a pair of bosons coupled to angular momentum 
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Similarly, the three-body interactions are given by 

L=0,2,3,4,6 i<j 
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The coefficients cl, and ^i^^ correspond to the 2 one-body, 7 two-body and 17 three-body matrix elements, 
respectively. They are chosen independently from a Gaussian distribution of random numbers with zero mean and 



variance v as 
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{Cu.Cl'^,^, ) = Sll' (1 + S,jyf )v^ , 

(eL'CL.,) ={eL<L„) = (a.,?L;,^.,) = 0. (6) 

The choice of the ensembles is such that they are invariant under orthogonal basis transformations. The variance of 
the Gaussian distribution sets the overall energy scale. The ensemble defined by Eq. (^) for Hk is called the /c-body 
random ensemble (fc-BRE) For two-body interactions H ^ H2 it reduces to the T(wo)BRE |l^Jl^. When the 
number of bosons is equal to the rank of the interactions N = k, the Hamiltonian matrix is entirely random and 
the ensemble coincides with the Gaussian orthogonal ensemble (GOE). For N > k the many-body matrix elements 
of iJfe are correlated via the appropriate reduction formulas and depend, in principle, on all random /c-body matrix 
elements. 



III. RESULTS 

In 1^ we used random one- and two-body interactions with TV = 16 to study the systematics of low-lying collective 
states in the IBM. Here we wish to study how these results depend on the boson number and the rank of the random 
interactions. 

We first analyze the dependence on the total number of bosons. Hereto we take the Hamiltonian H2 of Eq. (|^) 
with random two-body matrix elements. In all calculations we make 1000 runs. For each set of randomly generated 
two-body matrix elements we calculate the entire energy spectrum and the B{E2) values between the yrast states. In 
Fig. d we show the percentage of = 0+ ground states as a function of N (solid line). For N = 2 the Hamiltonian 
matrix is a real-symmetric random matrix. For each value of t he an gular momentum L the ensemble corresponds 
to GOE, whose level distribution is a semicircle with radius VMv^ and width \/{d + l)v^ [|l|. In this case, the 
percentage of ground states for a given value of L is determined by the dimension d of the Hamiltonian matrix: d = 2 
for L = 0, 2 and d — 1 ioi- L — 4. For 3 < < 16 the situation is completely different. The ensemble is now TBRE. 
The dominant angular momentum of the ground state is determined by the shapes of the level distributions as a 
function of the angular momentum, in particular by the tails, i.e. the higher moments, of the distributions (all have 
the same centroid). The distribution whose tail extends furthest is the most likely to provide the ground state. For 
a semicircular (GOE) or a Gaussian distribution (TBRE in the nuclear shell model [p2|-p^) the shape is completely 
detemined by the width. In these two cases, the dominance of = 0^ ground states can be correlated to the 
widths of the distributions However, for a system of interacting bosons the TBRE distribution of eigenvalues is 

neither semicircular (except for N = k) nor Gaussian [ p^ . There is no relation between the width (first moment) and 
higher moments of the distribution, which determine the dominant angular momentum of the ground state. Table || 
shows that the width increases with angular momentum, whereas the most likely value of the ground state angular 
momentum is = 0"*". In fact, the probability that the ground state has a certain value of the angular momentum 
is not really fixed by the full distribution of eigenvalues, but rather by that of the lowest one. Work is in progress to 
elucidate the form of these distributions in a schematic exactly solvable model . 

Despite the different shapes of the TBRE level distributions for fermions and bosons we find, just as in the fermion 
case [H-g], a dominance 60%) of ~ 0+ ground states in the IBM with 3 < TV < 16. This fraction is large 
compared to the percentage of — 0+ states in the model space (solid and dashed-dotted lines in Fig. p. The 
oscillations with maxima at = 3n (multiple of 3) are due to the 'unphysical' region of parameter space for which 
the energy ratio 

EiAt)-E{Ot) 
Ei2t)~E{0t)' 

is less than 1 (dashed line) corresponding to a level sequence 0^ , 4'^ , 2'^ , rather than 0^ , 2^ , i'^ for i? > 1 (dotted 
line) . The enhancement for N = 3n can be attributed to the existence of a 0+ state in which all d bosons are organized 
into riA = TV/3 triplets. This state has the U{5) quantum numbers \N,nd,v,nA, L >= |TV, TV, 0, TV/3, > and can 
become the ground state if the vibrational spectrum is turned 'upside down'. 

For the cases with a = 0"*" ground state we present in Fig. ^ the probability distribution P{R) of the energy 
ratio R of Eq. (^. This energy ratio has very characteristic values for the harmonic vibrator and the rotor, R — 2 
and R = 10/3, respectively. The Hamiltonian matrix of H2 depends on 7 independent random two-body matrix 
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elements. For small values of N there is little correlation among the matrix elements of H, and as a consequence 
the probability distribution P{R) shows little structure for iV = 3 (dashed-dotted curve). For increasing values of 
N there is a correspondingly higher correlation between the different matrix elements of H, which results in the 
development of two peaks in P{R). We first see the development of a maximum at i? ~ 1.9 for iV = 6 (dotted 
curve), followed by another one at i? ~ 3.3 for = 10 (dashed curve). For TV = 16 the probability distribution 
P{R) has two very pronounced peaks, one at R ^ 1.95 and a narrower one at i? ^ 3.35 (solid curve). These values 
correspond almost exactly to those for the harmonic vibrator and the rotor. The two maxima correspond to the 
two basic phases that characterize the collective region: a spherical one with i? ~ 2.0 and a axially deformed one 
with R ^ 3.3. There is no peak for 7-unstable nuclei (50(6) limit), since this requires that the matrix element 



of C2 



(fit X dt)(2) . (Jx d)(2) +h.C. 



vanishes identically, effectively corresponding to a zero-measure case for the 

random sample. Any other value of C212 (7^ 0) gives rise to an axially symmetric rotor p7| , p^ . 

In a second calculation we take the Hamiltonian H3 of Eq. (^ with random three-body interactions. Fig. |3| shows 
the same qualitative behavior as Fig. ^ although the peak structure is far less pronounced. For A'^ = 16 we see 
again two maxima at the vibrator and rotor values of the energy ratio R. The case of three-body interactions in the 
IBM is of special interest, since it can give rise to stable triaxial deformations (l9|, which are absent in the case of 
Hamiltonians including one- and two-body interactions only. We note, however, that in the neutron-proton version 
of the IBM, triaxial deformation can be obtained from Hamiltonians with one- and two-body interactions only [^ . 
Fig. H shows no clear sign of a 'triaxial' peak {e.g. a triaxially deformed rotor with 7 = 30° has R = 8/3), nor of a 
'7-unstable' one with R — 5/2. 



After these two model studies, we now turn to a more realistic case. It has been shown |10 that the phenomelogy 
of low-lying collective states in nuclei is well described by an IBM Hamiltonian consisting of both one- and two-body 
interactions 
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In order to remove the N dependence of the matrix elements of fc-body interactions, we have scaled Hk by ni=Ll^ + 
1 — i). In Fig. ^ we show the corresponding probability distribution P{R) of the energy ratio R of Eq. (Q) for 
different values of the number of bosons. The results are very similar to those of Fig. ^ which were obtained with 
pure two-body interactions. With increasing values of N the many-body matrix elements of H12 become increasingly 
correlated, which results in the development of two maxima in P{R). The curve for iV = 16 is identical to the 
calculation discussed in |^ . The occurrence of two basic phases for the collective region is further exemplified in Fig. ^ 
in which we plot the energy ratio R for the consistent-Q formulation of the IBM 
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with realistic values of the interactions, i.e. a positive d boson energy (e > 0) and an attractive quadrupole-quadrupole 
interaction (k > 0). The results in Fig. ^ are plotted as a function of the scaled parameters x = —2x/V7 and 
y ^ 4K(Af — 1)], which have been used as control parameters in a study of phase transitions in the IBM 

p8| , p2|j23| ) . For y = 1 we recover the vibrational or SU{5) limit of the IBM, whereas for y = and x = 1 one finds 
the rotational or SU{3) limit, and for y = and x = the 7-unstable or 50(6) limit. We clearly see two planes 
corresponding to i? 2.0 and R 3.3 respectively, which are separated by a sharp transitional region, in agreement 
with the observation in |2^,^ that the collective region is characterized by two phases (spherical and deformed) 
connected by a sharp phase transition. 

In order to investigate the effect of higher order interactions we now add three-body interactions to the Hamiltonian 
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In this case the Hamiltonian matrix depends on 26 independent random matrix elements (2 one-body, 7 two-body 
and 17 three-body). Therefore, for a fixed value of TV there is less correlation between the iV-body matrix elements of 
H123 than for H12, which results in broader peaks in the probability distribution P{R). A comparison of Figs. ^ and 
^ shows that the probability distribution P{R) behaves in a very similar way, and that the addition of three-body 
interactions does not change the results in a significant way. When N is sufficiently large compared to the maximum 
rank k of the interactions (2 and 3, respectively) the results become independent of k. 

This result is qualitatively very similar to that of [ p^ , in which the transition from a Gaussian to a semicircular 
level distribution was studied in the nuclear shell model for a fixed particle number N — 7 with increasing values of 
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the rank 2 < k < 7. The characteristic features of the ensemble depend on the ratio of the number of particles and 
the rank of the interactions. For N sufficiently large compared to k there is a saturation, and the properties of the 
ensemble no longer depend on k. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied global properties of low-lying collective levels using the interacting boson model with 
random interations. In particular, we addressed the dependence of the dominance of = 0+ ground states and the 
occurrence of vibrational and rotational band structures on the boson number N and the rank k of the interactions. 

Just as for the nuclear shell model it was found that despite the randomness of the interactions (both in size and 
sign) the ground state has — 0+ in approximately 60 % of the cases. The oscillation in the percentage of = 0+ 
ground states with N was shown to be entirely due to cases in which the level sequence is given by 0+, 4+, 2+ (i? < 1). 
For the cases with R> 1 there is a very smooth dependence on TV. 

The vibrational and rotational band structures appear gradually as N/k increases. For N ^ k there is little 
or no evidence for such bands. As N grows we first see evidence for the development of vibrational structure, 
followed later by the appearance of rotational bands. If N increases further these band structures become more and 
more pronounced. Essentially the same behavior is found for random two- and three-body interactions. In realistic 
applications to collective nuclei the IBM Hamiltonian consists of a combination of one- and two-body interactions. A 
study with random ensembles of one- and two-body interactions shows similar results to the case of pure two-body 
terms. The inclusion of random three-body interactions does not significantly change the basic features. 

In conclusion, we find that the dominance of — 0^ ground states and the occurrence of vibrational and rotational 
features are independent of the boson number, as long as A'^ is sufficiently large compared with the maximum rank 
of the interactions. We can conclude that these features represent general and robust properties of the interacting 
boson model space, and are a consequence of the many-body dynamics, which enters via the reduction formulas for 
the A'^-body matrix elements of fc-body interactions (angular momentum coupling, coefficients of fractional parentage, 
etc.). Since the structure of the model space is completely determined by the corresponding degrees of freedom, 
these results emphasize the importance of the selection of the relevant degrees of freedom. In this context, a relevant 
question is whether vibrational and rotational collective behavior can be directly observed in the shell model with 
random interactions if an appropriate truncation of the (shell model) Hilbert space is carried out. 

It is important to stress that these properties do not arise as an artefact of a particular model of nuclear structure. 
In empirical studies of the low-lying collective states of medium and heavy even-even nuclei very regular and robust 
features have been observed, such as the tripartite classification into seniority, anharmonic vibrator and rotor regimes 
and the systematics of excitation energy and Ml strength of the scissors mode [ p4[ . 

Finally, we remark that the use of random interactions to study the generic behavior of low-lying states has also 
found useful applications in many-body quantum systems of a different nature, such a quantum dots or small metallic 
particles [ p5[ . 
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TABLE I. Percentage of ground states with angular momentum L and corresponding widths of TBRE level distributions. 
The results are obtained for 1000 runs and iV = 16 bosons. The widths are divided by N{N — 1). 
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FIG. 1. Percentage of = 0"'' ground states as a function of the boson number A'' for H — H2 for which the energy ratio 
of Eq. (0) is given by < ii < 1 (dashed line), i? > 1 (dotted hne) and R> (sohd line). The dashed-dotted line shows the 
percentage of — 0"*" states in the model space. 
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FIG. 2. Probability distributions P{R) of the energy ratio R of Eq. ^ with J P{R)dR = 1 in the IBM with random two-body 
interactions for TV = 3 (dashed-dotted), 6 (dotted), 10 (dashed) and 16 (solid). 
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FIG. 5. The energy ratio 7? of Eq. (M) as a function of x and y in the consistent Q-formulation of the IBM. 




